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Abstract: It has been conjectured that the classical dynamics of M theory is equiv- 
alent to a null geodesic motion in the infinite-dimensional coset space Eio/K(Eio), 
where K(Eio) is the maximal compact subgroup of the hyperbolic Kac-Moody group 
Eio- We here provide further evidence for this conjecture by showing that the lead- 
ing higher order corrections, quartic in the curvature and related three-form dependent 
terms, correspond to negative imaginary roots of Eio- The conjecture entails certain 
predictions for which higher order corrections are allowed: in particular corrections of 
type R M (DF) N are compatible with E w only for M + N = 3k + 1. Furthermore, the 
leading parts of the R 4 , R 7 , • ■ ■ terms are predicted to be associated with singlets under 
the alio decomposition of -Eio. Although singlets are extremely rare among the alto- 
gether 4 400 752 653 representations of slio appearing in Eio up to level i < 28, there 
are indeed singlets at levels i = 10 and £ = 20 which do match with the R 4 and the 
expected R 7 corrections. Our analysis indicates a far more complicated behavior of the 
theory near the cosmological singularity than suggested by the standard homogeneous 
ansatze. 



1 Introduction 

The analysis a la Belinskii, Khalatnikov, Lifshitz (BKL) [ 1 1 of generic cosmolog- 
ical solutions of D = 11 supergravity [2] in the vicinity of a spacelike singularity 
has revealed a connection with billiard motion in the fundamental Weyl chamber 
of Eio (implying chaotic oscillations of the metric near the singularity) QJHI|5)- 



We recall that E\q is a rank- 10 infinite-dimensional hyperbolic Kac-Moody alge- 
bra 1 whose root lattice is the canonical hyperbolic extension of the root lat- 
tice of the largest exceptional finite-dimensional Lie algebra E%, and the unique 
even self-dual Lorentzian lattice Hi 9 [7]. The cosmological billiard describing the 
asymptotic behaviour near a spacelike singularity is based on the identification of 
the ten diagonal metric degrees of freedom of D = 11 supergravity with the ten 
non-compact directions of a Cartan subalgebra (CSA) of E\q. This inuiguing link 
between D = 11 supergravity and E\q was deepened in [8 9] where it was shown 
that the bosonic equations of motion of D = 11 supergravity at some given spatial 
point, when restricted to zeroth and first order spatial gradients in the metric and 
the three-form, can be matched with the equations of motion of a one-dimensional 
E\o/ K{E\q) (T-model restricted to levels £ < 3. In terms of the heights of the 
roots of £10 involved in these correspondences, the Weyl-chamber billiard retains 
only roots of height one (i.e. simple roots), while the work of [ 8 9 1 has established 
the correspondence with (null) geodesic motion on the coset space E\q/ K{E\q) 
up to height 29 included. These results underline the potential importance of E\q, 
whose appearance in the reduction of D = 11 supergravity to one dimension had 
been conjectured already long ago in fTOllTTl . as a candidate symmetry underlying 
M theory. A similar, but conceptually different, proposal was made in fT2l . where 
En (or some even larger symmetry containing E\\ lTP3l ) has been suggested as a 
fundamental symmetry of M theory. Let us also note that links between the dynam- 
ics of gravitational theories and geodesic motions on certain Lorentzian signature 
spaces were noticed long ago lfl4l[T5l . For recent recent work in this vein in the 
context of homogeneous cosmologies (but not related to E\q), see fT6l[T7l . 

The substantial increase in the height of the roots entering the correspondence 
between the two dynamics led to the key conjecture of |8), according to which 
the one-dimensional bosonic E\q/ K(Eiq) cr-model is equivalent (or 'dual') to the 
bosonic sector of D = 11 supergravity, possibly augmented by further M-theoretic 
degrees of freedom. To establish this conjecture, one faces two challenges, namely 
to understand 

• how the higher order spatial gradients are realized in the cr-model, and how 
the full (untruncated) bosonic sector of D = 11 supergravity is thereby em- 
bedded into the E\q/K{E\q) coset space dynamics; and 

• the meaning and significance of the imaginary (lightlike and timelike) roots 
of E\q and their relation to M theory degrees of freedom and M theory cor- 
rections of the 'low energy' D = 11 supergravity action. 

'For simplicity of notation, E^q denotes both the group, and its associated Lie algebra. 
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Concerning the first challenge, it was shown in 1 8 1 that Eiq contains three 
infinite towers of Lie algebra elements, which possess the correct structure for rep- 
resenting the higher spatial gradients of the metric, and of the 'electric' and 'mag- 
netic' components of the four-form field strength. However, it is still unknown 
how this identification could work in detail. The progress reported in the present 
paper concerns the second challenge. Namely, we shall show that certain (partially) 
known higher-order M theory corrections to the two-derivative, 'effective' D = 11 
supergravity Lagrangian can be related, in a certain approximation which we shall 
explain, to special imaginary roots of E\q. Our 'botanical' approach in search- 
ing for new connections between M theory and the hyperbolic Kac -Moody algebra 
.Eio is based on the working hypothesis that the Eiq /K (Eiq) a-model of ^ not 
only describes M theory degrees of freedom beyond those of D = 11 supergravity, 
but in addition contains hidden information about M theory corrections to the low 
energy effective theory (i.e. D = 11 supergravity) at arbitrarily high orders. More 
specifically, we shall consider the leading (8th order in derivatives) corrections to 
the usual supergravity Lagrangian and exhibit their connection with certain imag- 
inary roots of E\q. Different routes for interpreting the imaginary roots of E\q in 
terms of brane dynamics have been suggested in [ 18 1, and, in the framework of the 
E n proposal of in 111911201 . 

In [ 5 1 it was shown that the gravitational Hamiltonian, or more precisely the 
Hamiltonian constraint at a given spatial point, can be written in the form 



H(P a ,ir a ,Q,P) = N 



-G ab 7r a 7r b + Y,ca(Q,P, d(3, 8Q) exp ( - 2w A (f3)) 
1 A 

(1) 

with the rescaled lapse N = N/ yjg, where g is the determinant of the spatial met- 
ric. Here (3 a , 7r a are the canonical variables corresponding to the diagonal (spatial) 
metric degrees of freedom, and G ah the (Lorentzian) 'superspace' metric induced 
by the Einstein-Hilbert action. (Q, P) denote the remaining canonical degrees of 
freedom associated to off-diagonal metric and various matter degrees of freedom, 
and dj3, dQ their spatial gradients. The sum runs over various exponential po- 
tential 'walls', which are indexed by A. The above form of the Hamiltonian is 
in particular valid for the bosonic sector of D = 11 supergravity, in which case 
0,6 = 1,..., 10. 

On the other hand, for the one-dimensional 'geodesic' <j-model over the infinite 
dimensional coset space E\q/ K(E\q), an analogous expression for the Hamilto- 
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nian constraint was also derived in [5 ], viz. 



p) = n ^G a \ a 7T b + Yl E (n a ,s(^,p)) 2 exp(-2a(/3)) 



mult (a) 



H(p a ,ir a ,v, 



a€A + s=l 



(2) 



where (/3 a , 7r a ) now denote the ten diagonal CSA degrees of freedom of E\q, and 
{v,p) stand for infinitely many off-diagonal (Iwasawa-type) canonical variables, 
on which the quantities Il QjS depend. The metric G ab now denotes the restriction 
of the unique invariant bilinear form on E\q to the CSA J6j. This metric happens 
to be identical with the metric appearing in for the bosonic sector of D = 1 1 
supergravity; its explicit form is given in d26l below. The sum on the r.h.s. of © 
ranges over all positive roots a of E\q with their multiplicities [= mult(a)], and 
the number a(/3) is the result of applying the linear form (=root) a to the CSA 
element 0. The formula © is actually valid for any G/K(G) cr-model over an 
indefinite Kac-Moody group G with 'maximal compact subgroup' K{G). Varying 
© w.r.t. the (new) lapse n implies that the resulting solution is a null geodesic on 
this infinite dimensional coset manifold. 

The similarity of Q and © is obvious, but the coincidence of the two expres- 
sions and of the associated bosonic equations of motion has so far only been es- 
tablished for a limited number of terms. The first check consists in matching those 
terms in (0 and © which dominate the dynamics near a cosmological singular- 
ity. In this leading BKL approximation one can prove [5| that the 'off-diagonal' 
degrees of freedom (Q, P) in Q and (u, p) in (|2j both 'freeze' as one approaches 
the singular initial hypersurface T = 0. [Here T denotes the proper time. We shall 
reserve the letter t to denote the coordinate time in the gauge N = 1, i.e. N = ^fg. 
This coordinate time goes to +oo as T — > 0.] The precise form of the coefficient 
functions ca{Q, P,d[5,dQ), which is very complicated, does not matter in this 
BKL limit. What matters is that the coefficients of the leading contributions are 
all non-negative. That is, we have ca' > 0, where the primed index A' labels the 
leading terms. This implies that, for large (3, the Hamiltonian Q takes the limiting 
form 



where the sum is only over the 'dominant walls'. It was also shown, that in this 
limit the exponential (Toda) walls appearing in © can be replaced by infinitely 
high ('sharp') walls. Technically, this means that one can replace the exponential 
functions in (|3} by an infinite step function © [5]. The billiard then takes place 
in the chamber defined by 0(— 2wa'{/3)) = 0, or equivalently, wa'(P) > 0. An 
analogous argument applies to the coset dynamics ©. The freezing of all the 



Woo(/3 a ,7r a ) = -G ab 7r a vr 6 + ^c A ,exp(-2^(/3)) 



(3) 



.4' 
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off-diagonal degrees of freedom u,p implies that the combinations Ii a ^ s {v,p) also 
freeze near the singularity, so that the Hamiltonian © takes the limiting form (in 
the gauge n = 1) 

#oo(/? a , vr a ) = ^G ab Tr a ir b + < ex P ( " ^O 9 )) W 

where the sum over the positive roots of E\o can be asymptotically restricted to 
the subset of leading positive roots a{. It is easily seen that these 'leading' positive 
roots are nothing but the simple roots of Eiq. 

The equivalence between the two limiting dynamics © and © is now a con- 
sequence of the fact that the set of dominant supergravity walls wa> (/?) coincides 
with the set of simple roots on{0) of E%q E). The validity of this equivalence was 
extended in [8 9| from the height one simple roots to height 29 in E\q roots . 
At this order, one keeps more exponential terms in both Q and ©. Now, the 
precise form of the supergravity coefficient functions ca(Q, P, . . .) does matter, 
and was shown in |8] |9) to correspond precisely to the form of the coset functions 
n, aiS (v,p), with a 'dictionary' relating the variables of both dynamics. This match- 
ing between the two Hamiltonians was checked at the level of the equations of mo- 
tion. It is still unclear whether (and how) one can extend the dictionary relating the 
dynamical variables in both models so as to prove the equivalence between the two 
Hamiltonians © and © beyond height 29. It is quite possible that this dictionary 
becomes spatially non-local in the gravity variables, and that the simple identifica- 
tion assumed above between the diagonal gravitational degrees of freedom /3 a , 7r a 
and the Eiq CSA variables entering $2% must be modified beyond height 29. Such 
non-trivial changes of variables may be needed, for instance, in order to reconcile 
the fact that the only negative contribution in (0 comes from the CSA (the first 
term on the r.h.s.), whereas in the expansion of the gravity Hamiltonian © there 
may arise negative contributions also from (subleading) terms with ca < — 
as follows from considering, for instance, homogeneous spatial geometries with 
curvature of arbitrary sign. 

Our aim here is to investigate what types of wall forms can be formally asso- 
ciated with higher order terms (linked to quantum corrections in M-theory), no- 
tably the leading ones of the type R 4 ,R 2 (DF) 2 , R(DF) 3 and (DF) 4 , as well 
as the topological Chern Simons term A3ARARARAR. Remarkably, we 
shall find that these terms correspond, in leading approximation, to time-like imag- 
inary roots of E\q. Recall that the 'length' of a root is calculated by means of 
a quadratic form G ab , which is the same as the one appearing in (0 (and ©), 

2 The height of a positive root a — n i a i ( n i > 0) is defined as ht(ai) = Yl n i I- 1 1 [6). 
Correspondingly, the height of a negative root a — — Y, THOti (rii > 0) is defined as — Yl n i- 



5 



and which, for indefinite Kac-Moody algebras, is Lorentzian. Spacelike roots, for 
which G ab a a at, = a 2 > 0, are called 'real roots' in the mathematical literature, 
while light-like or time-like ones, for which a 2 < 0, are called 'imaginary roots' 
0. All the roots up to height 29, that have been explicitly checked to match be- 
tween the two Hamiltonians Q and © turned out to be real and positive. By 
contrast, we will find here that the leading roots associated with 8th-order deriva- 
tive corrections are, at once, imaginary {a 2 < 0) and negative ( i.e. of the form 
a = — nidi with m > 0). For the leading R 4 corrections, they lie deep inside 
the root-space lightcone, with a 2 = —10, at height ht(a) = —115, and at level 
1 = -10. 

The appearance of exponential walls exp ( — 2w{(3)) associated with negative 
roots in a version of the gravitational Hamiltonian Q with higher order corrections 
included, might seem to be in plain contradiction with the structure © of the coset 
Hamiltonian, which contains a sum restricted to positive roots a S A + . However, 
one must remember that the explicit form of the coset Hamiltonian strongly de- 
pends on the parametrisation used for a general coset element V G E\q/K{E\q). 
The form © results from an Iwasawa parametrization V = VB, where V is 'di- 
agonal' {i.e. in the Cartan torus of E\q) and B of Borel, or 'upper triangular' 
type, i.e. obtained by exponentiating a sum of positive root generators. If instead 
we were to employ a more complicated parametrization of a generic coset element 
V £ E\q/ K{E\q), namely involving, besides the exponentiation of infinitely many 
positive root generators, the exponentiation of some 3 negative root generators, one 
would end up with an Hamiltonian of a more general form than The latter 
would involve some mixture between positive-root walls exp ( — 2aA(/3)) (with 
a a € A + ), and negative-root ones exp ( + 2ob(/3)) (with as G A+). Actually, 
the mixtures of positive-root walls and negative-root ones arising in such non-Borel 
parametrizations are rather complicated, and do not reduce to a decoupled sum of 
normal exponentials exp ( — 2a^(/3)), and inverse ones exp ( + 2ob (/?)). A non- 
Borel parametrization for the root a>B tends to generate superpositions of effects 
proportional to both exp ( + 2ob(/3)) and exp ( — 2ob (/?))• Therefore, the wrong- 
sign walls ~ exp ( + 2ob(/?)) can be expected to dominate over the usual-sign 
ones only in some special corner of the dynamical space. Such an approximate 
dominance is enough for our purposes here, because we shall see below that the 
negative root walls associated to R 4 terms can only be meaningfully considered in 
some intermediate domain of the dynamical evolution. 

No doubt there is still a long way towards the ultimate goal of proving the 

3 We note that such 'mixed' non-triangular parametrizations of infinite dimensional Kac-Moody 
coset spaces have so far not been considered in the literature. It seems likely that, for the exponen- 
tiation to be well defined for hyperbolic Kac-Moody groups, one can only have a sparse number of 
negative root generators mixed with an infinite number of positive root ones. 
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equivalence between the (essentially unique) E\q coset dynamics © and the full 
effective (bosonic) supergravity action (involving the quantum effect of further M- 
theoretic degrees of freedom), but we hope that the results reported here will allow 
one to make new and stringent tests of our basic conjecture, and to make new pre- 
dictions concerning both the structure of higher-order corrections in supergravity 
and M-theory, which are completely inaccessible by conventional methods, and the 
existence of special slio-representations in the level decomposition of E\q. 

2 Asymptotics of frames and scale factors 

When working with D = 11 supergravity fields, we use Lorentz (tangent space) 
tensors throughout, with the 'mostly plus' metric r\ AB = ( — !-•••+), and the 
following index conventions 

Flat spacetime indices : A, B, C, ... € {0, 1, ... , 10} 

Flat spatial indices : a, b, c, ... G {1, ... , 10} (5) 

When needed, the coordinate spacetime indices will be denoted as M, N, P, Q and 
the coordinate space ones as m, n, p, q. 

The Lorentz covariant derivative is defined in the usual manner on any Lorentz 
vector Va via 

D A V B :=d A VB+0JAB C Vc (6) 

where 8a is a (non-commuting) frame derivative linked to the usual (commuting) 
spacetime derivative 8m through 8a = Ea M 8m, where Ea M is the elfbein. The 
spin connection is given by the standard formula 

t-^ABC = 2 {^ABC - ^BCA + ^CAB^j = ~^A CB (7) 
in terms of the coefficients of anholonomicity 

n A B C := Ea M E b n {8mE n c - 8 N E M C ) = -Vba C (8) 

where Em is the co-frame, i.e. the inverse of the elfbein Ea ' ■ Expressed in flat 
indices, the Riemann tensor is 

Rabcd = 9aujbcd — Qbloacd + ^ab E ^ecd 

+^ac e ^bed — ^>bc e ^>aed (9) 

Likewise, for the 4-form field strength we will also use flat indices such that 

Fabcd = Ea m E B N E c p E D ® F MNP q , etc. 
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For the further analysis, we will rely on the following (zero-shift) 1+10 split of 
the elfbein 

with the spatial (inverse) zehnbein e rn a , and adopt the gauge 

N = ^ = dete m a (11) 

As explained in 0, and as is also obvious from the definition of N in Q, one 
consequence of this choice is that the leading kinetic term in © simplifies to oc 
G ab ir a ir b . Splitting indices into A = and A = a, where denotes a 'flat' (proper) 
time index and a a flat spatial index, the coefficients of anholonomicity become 

^abc — [a 6ty d m €. nc , 

tt obc = N~ l e b n d t e nc , 

= ^ooa = -e a m N~ l d m N = -e a m e- l d m e (12) 

with all other coefficients of anholonomicity vanishing. From the above formulas 
we obtain in particular 

w fec = N~ l e [b n d t e nc } , u a w = N~ l e {a n d t e nb) (13) 

where (...) and [...] denote symmmetrization and antisymmetrization with strength 
one, respectively. 

In the sequel we will analyze the leading asymptotic behaviour of various fields 
near T = as functions of the logarithmic scale factors (3. To this aim we intro- 
duce, as an intermediate object, the Iwasawa co-frame 9 m a = A/" m , where A/" m is 
the upper triangular matrix entering the Iwasawa decomposition of the metric g mn 
that was introduced in 0. [This triangular matrix has ones on the diagonal, so 
det = 1.] In terms of this intermediate frame, the ten-dimensional spatial Cartan 
moving frame (zehnbein) is given by (no summation on a) 

e m a = e-? a 6 m a (14) 

in each spatial section t = const, (recall that we use a zero-shift gauge). The line 
element then reads 

ds 2 = -(Ndt) 2 + Y,(e-^n 2 d5) 

a 

The decomposition (THT i of the zehnbein is most useful in analyzing the small T be- 
havior of the gravity model. More specifically, the dynamical variables (3 a and 9 m a 
behave very differently as T — > 0. While the frame 8 m a was shown to 'freeze' (i.e. 
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to have a well-defined limit) as T — > 0, the logarithmic scale factors f3 a (whose 
time-derivatives are proportional to the Kasner exponents) generically diverge to- 
ward +00 near the singular initial ('big bang') spatial hypersurface T = 0, while 
exhibiting chaotic oscillations of BKL type. [Note, in passing, that we never as- 
sume the /3's to be only functions of time t. All our supergravity dynamical vari- 
ables are functions of all eleven coordinates (t,x m )]. With these notations, the 
gauge condition (fTTl implies 

N = exp f-J^J d ete = e~ a (16) 

while the D = 11 volume density behaves as 

E = N^ = e- 2a (17) 

Here a denotes the sum of all logarithmic scale factors, 

10 

<r:=5> a (18) 

0=1 

This combination will play an important role in the remainder. As T — > 0, the 
volume factor *Jg = e~ a tends toward zero, and a tends to plus infinity. Similarly, 
one can exhibit the behavior of various other fields near the singular initial hyper- 
surface by factoring them into a part that 'freezes' as T — > and is independent 
of the scale factors, and another, dependent on the /3's, which diverges and may 
oscillate near T = 0. 

Making use of (fT4l we can compute the asymptotic behavior of the coefficients 
of anholonomicity and the spin connection. Suppressing inessential prefactors, we 
find for the connection components with one (proper) time index 

wok ~ 

ooabo = -N-H ab d t p a + ... = -5 ab e +a d t p a + ... (19) 

where the symbol ~ indicates that we keep only the leading contributions (from 
(fl2l it is easily seen that that the components with two time indices are sub- 
dominant). For the purely spatial components, we get for the anholonomicity co- 
efficients 4 

Oq be = e^+^ fW + Sa^'B^ - 5 bc e^B a p c (20) 

4 The spatial indices in these and other formulas below always refer to the basic orthonormal 
frame e m a , whereas the barred quantities associated to the intermediate (non-orthonormal) frame 
9m a are only used as auxiliary objects. 
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where d a = a m d m , and 

^abc '■= '26[a m Qb\ n dm8nc (21) 

Hence the spatial components of the spin connection are given by a sum of terms: 

+ 5 ac e pb d b (3 c - 5 ab e? a d c p c (22) 

It is important that the barred quantities 0, abc and the frame derivatives B a all pos- 
sess finite limits as T — » 0; they correspond to the 'frozen' components in the 
language of [8]. Hence, for T ~ 

io a b c oc +/3 ~P + permutations of the indices a, b, c (23) 

where we may have a = b or a = c (but b ^ c as to abc = —^>acb)- We will make 
use of these formulas in section 5, when determining the asymptotic behavior of 
the curvature and other 'composite' quantities of interest. 



3 Dominant walls and E m roots 

In this section, we recall some basic facts and useful formulas concerning the Lie 
algebra E\q and its roots. A crucial fact here is the identification of the logarithmic 
scale factors appearing in (fUT l with the CSA degrees of freedom of E w . This will 
enable us to re-interpret the exponential walls exp[— 2w(f5)\ in the gravitational 
Hamiltonian in terms of the exponential (Toda) walls exp[— 2a(/3)] entering the 
E\q cr-model ®. As mentioned already, the scalar product on the space of scale 
factors induced by the Einstein-Hilbert action (for D = 11) coincides with the 
restriction of the invariant bilinear form on E'io to its CSA [5|. This scalar product 
is given by 

10 / 10 \ 2 

(f3\P) = G ab p a (3 b = ^ (/T) 2 - hT/H (24) 

a=l \a=l J 

It is indefinite (Lorentzian) due to the lower unboundedness of the Einstein-Hilbert 
action under variations of the conformal factor. As generally shown in 0, in the 
limit T — > the dynamics of the scale factors is well approximated by a billiard. 
The region defining this billiard is contained in the forward lightcone in /?-space, 
and bounded by 'sharp walls'; the latter are hyperplanes defined as the zeros of 
certain linear forms ('wall forms') 

10 

= J>a/? a (25) 

a=l 
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The metric on the dual space of wall forms is the contravariant form of d24l 

10 i / 10 \ 2 

(p\ P ) = G ab PaPb = E - q < 26) 

o=l \a=l / 

and this is also the metric appearing in (0 and ©. Below, we will also use the 
notation 

w = (pi,...,p 10 ) (27) 

to designate a particular wall form d25l . In [ 5 ] we have explained how to determine 
these wall forms from the matter coupled Einstein-Hilbert action in a canonical 
formulation (see also the following section). A main result was that for D = 11 
supergravity, each wall form w can be associated with a particular real root of E\o, 
in accordance with the fact that the root space is dual to the CSA, which itself is 
identified with the space of scale factors f3. In particular, the dominant walls are 
associated with the simple roots. In the above basis d27l . the simple roots of E\q 
are 



ao = 


(1,1,1,0,0,0,0,0,0,0) 


Ot\ = 


(-1,1,0,0,0,0,0,0,0,0) 


(*2 = 


(0,-1,1,0,0,0,0,0,0,0) 


a 3 = 


(0,0,-1,1,0,0,0,0,0,0) 


£*4 = 


(0,0,0,-1,1,0,0,0,0,0) 


Q 5 = 


(0,0,0,0,-1,1,0,0,0,0) 


«6 = 


(0,0,0,0,0,-1,1,0,0,0) 


a-j = 


(0,0,0,0,0,0,-1,1,0,0) 


a 8 = 


(0,0,0,0,0,0,0,-1,1,0) 


Q 9 = 


(0,0,0,0,0,0,0,0,-1,1) 



and one easily checks that they generate the Eiq Dynkin diagram via the scalar 
product d26l i. In the billiard picture, the dominant wall forms are therefore 

a (P) = /3 1 + /3 2 + /? 3 
aiO0) = /3 2 -/? 1 

ag(0) = P w - (3 9 (29) 

The first simple root corresponds to a leading 'electric wall', while the remain- 
ing nine simple roots correspond to the leading 'symmetry walls'. All other walls 



w ao (/3) = 

Wai(P) = 

w ag (f3) 
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(electric and symmetry, as well as the magnetic and gravitational walls) are sub- 
dominant and can be written as positive linear combinations of these leading walls. 
Consequently, the billiard takes place in the fundamental Weyl chamber, namely 
the 'wedge' defined by the conditions oti((3) > 0, or 

(3 1 + (5 2 + (3 3 > , (3 1 < I3 2 < . . . < (3 10 (30) 

Every positive root of Exo can be expressed as a linear combination of simple roots: 

9 

a = £ag + m? oij = [I ; m 1 , m 2 , m 3 , m 4 , m 5 ,mP,m 7 , m 8 , m 9 ] (31) 
i=i 

where, of course, £, m? > 0, and the height of a is ht (a) = £ + Ylj mJ • The 
brackets ([• • •]) denote the components of a root in the basis of simple roots, and 
are to be distinguished from the use of round parentheses used to denote the com- 
ponents of a root in the 'coordinate basis' of the /?'s. The integer £ is referred to as 
the 'Ag level', or simply the 'level' of the root a. (Different 'levels' can be defined 
w.r.t. other subalgebras of Exq such as Dg fZ2\ or Eg [23]). This level provides a 
grading of the Kac-Moody algebra Exo, allowing a decomposition of the Exo Lie 
algebra into an infinite tower of slxo representations, which has been worked out 
up to level £ < 28 in l24l . The two bases d27t and (ITU are related by 

px = £ — m 1 

p2 = £ + m 1 — m 2 

P3 = £ + m 2 — m 3 

Pi = to 3 — m 4 , ... , pg = m 8 — m 9 , pxo = m 9 (32) 

Roots corresponding to wall forms associated with particular higher order terms 
will appear at a fixed level £. For the associated roots we will usually give the one 
corresponding to the root of lowest height within a given slxo multiplet, for which 
the integers p a are ordered according to 

Pi > P2 > ■ ■ ■ > Pio (33) 

Because this is also the lowest weight of the corresponding SL (10) representation, 
we can read off the associated Dynkin labels directly from the tables of l24l . as 
well as the outer multiplicity fj,, i.e. the number of times this representation occurs 
in the decomposition of Exo- Starting from the lowest weight, other weights in 
the representation can be reached by successive addition of the 'symmetry roots' 
ax, ■ ■ ■ , ag. From d28l it is evident that this operation simply amounts to a permu- 
tation of the p^s (and in in terms of the scale factors (3 a to a replacement of the 
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spatial indices a by other spatial indices). The highest weight, which corresponds 
to a root of maximal height, is reached when the order above has been inverted to 
Pi < ■ ■ ■ < Pio- 

As an example let us briefly recall how the gravitational billiard walls arise 
by expanding out the Hamiltonian constraint following from the Einstein-Hilbert 
action, as described in detail in [5 1, and how they can be identified with particular 
E\q roots. Using d20b . they are obtained by considering the following contribution 
(no summation on repeated indices) 

En abc n abc = e~ 2(7+2l3a+2l3b -^ c n abc n abc + subleading terms (34) 

Because the walls always appear via exponentials e~ 2w ^\ the relevant linear form 
is 

w 9™ v = a-(3 a -(3 b + (3 c , with a ^ b (35) 

Choosing the indices (a, b, c) = (9, 10, 1) so as to obtain the root of lowest height 
in this multiplet [with Dynkin label (100000010)] we get 

<£i = (2, 1,1, 1,1, 1,1, 1,0,0) ^ a = [3; 1,3,5,4,3,2,1,0,0] (36) 

In particular the associated root is real {a 2 = +2) and positive in our conventions. 
Consequently, for D = 11 supergravity, the above gravitational walls arise at level 
I = 3 and are thus 'behind' the level I = 1 electric walls defined by the simple root 
ao- This contrasts with pure gravity, where the leading gravitational wall always 
corresponds to a simple root of the algebra AE n j4). 

4 Higher order corrections: what is known 

Still as a preparation, we here summarize briefly what is known about the 8th order 
derivative corrections to the eleven-dimensional supergravity action, which are of 
the type R'\ R 2 (DF) 2 , R(DF) 3 , (DF) 4 , • • • , F 8 , where R stands for the cur- 
vature (Riemann) tensor, and F = dA is the four-form field strength (there are 
no R 3 DF terms because any possible contraction of the indices will involve an 
e-tensor, and vanish by the Bianchi identity for either R or DF). In addition, 
there are separate 8th-order derivative terms of Chern Simons type proportional to 
A3R 4 required for the cancellation of anomalies [25|. The supersymmetric com- 
pletion of these 8th-order terms was studied in [26 1 and, in a superspace formu- 
lation, in [27 28 1 (see also the review [29 1). The presence of such corrections 
can be inferred (i) from (string) one-loop corrections to D = 10 type IIA string 
amplitudes ll30l OTl . which are similar in form to (a') 3 corrections at string tree 
level 1321 l33l l34l 1351 ; (ii) from considering M-theory corrections at one loop 



13 



l36l l37l l38l l39l . (Hi) from the computation of 4-graviton amplitudes in super- 
membrane theory by use of supermembrane vertices |4"U1I4T1I4"2"1: or (iv) from the 
structure of divergences in D = 11 supergravity 1431 . However, none of these 
approaches gives the complete 'reduced' {i.e. maximally simplified) kinematical 
structure of all the (bosonic) 8th-order derivative corrections. The paper [43] does 
contain explicit expressions for all terms oc R M (DF) N (M + N = 4), but in an 
'unreduced' form involving the Bel-Robinson tensor, which for our present pur- 
poses is rather cumbersome 5 . Moreover, none of these papers gives expressions 
for the terms involving the undifferentiated 4-form F, such as F s . 

We note that the D = 11 results should be compatible with those of the IIA 
theory, whenever the D = 10 results can be 'lifted' to D = 11. In this case 
one must worry about the possible existence of terms that algebraically vanish in 
D = 10, but not in D = 11. As the number of independent invariants made out 
of the Weyl tensor is = 7 when D > 8 1441 . this is not a concern for the curvature 
terms. However, the situation is more subtle for invariants built from the covariant 
derivative DF of the 4-form field strength, whose number does depend on the 
dimension. For instance, the number of (DF) 4 invariants is found to be = 23 
for D = 9, and = 24 for D = 11, as can be conveniently checked by use of the 
computer algebra package LiE [45 1. For D = 10, we would have 29 invariants with 
a 4-form field strength; however, there are only 15 invariants that can be built from 
the NSNS 3-form field in IIA supergravity (which means that we cannot simply lift 
the results of [35] to eleven dimensions). 

Yet another subtlety one must keep in mind is that the lightcone results do not 
necessarily agree with results of a 'covariant' calculation. A case in point is the 
'Euler-Lovelock' term which we will encounter below: this term does not show up 
in the lightcone calculation of the 4-graviton amplitude, hence is not determined, 
whereas it does appear in the covariant approach of l39ll26l . However, this term 
should start to contribute to on-shell amplitudes from 5-point amplitudes onwards, 
whence it should become visible in a lightcone computation of the 5-graviton am- 
plitude. 

We are here only interested in on-shell counterterms. That is, we shall ignore 
all terms proportional to the equations of motion because they can be absorbed into 
redefinitions of the basic fields. For this reason, the Riemann tensor Rabcd can 
be replaced by the Weyl tensor Cabcd, since all terms containing the Ricci tensor 
or Ricci scalar can be absorbed (modulo F-dependent terms) into a redefinition 
of the metric by use of the equations of motion. For the DF type terms, we can 



5 Beware of a missing factor 12 in the third term, oc P 2 , on the r.h.s. of Eq. (5a) of the first 
reference in 1 4-3 (J. Plefka, private communication). 
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similarly discard the trace and the fully antisymmetric combinations because 



D Fabcd ~ 



D[aFbCDE] 







(37) 



by the equations of motion and the Bianchi identity (we neglect terms oc F 2 on the 
r.h.s. of the equation of motion; we shall see later that they are indeed subdomi- 
nant). Thus we must only deal with the irreducible representations of the rotation 
group represented by the following Young tableaux: 



Cabcd 



DaFbcde 



(38) 



We first focus on the terms quartic in R (or C), because, as we shall see, they 
are the dominant ones within our analysis. From Refs. l39l |26 ] these terms are 
given by the expression (modulo a positive coefficient) 



£ (4)tot = E ( Jq _ 2l2) 

where E = det E M A = N^fg, 



4 ' 1 E 1 - V \ Z 



Jo :— tgtgR + -Es = X 



(39) 



(40) 



and 



To 



1 



-E 8 +2e n A 3 



tiR 4 



-(tvR 2 ) 2 



= --Z+2e u A 3 



tvR 4 



(tvR 



2^2 



(41) 

Here, we use the same condensed, index-free notation as in [39], while the notation 
X = t 8 t$R 4 and Z = -2E S is used in (26). 6 

In particular, tg denotes the 8-index tensor that enters the kinematic factor of 
both tree-level and one-loop 4-particle string amplitudes 1301 l46l . The normaliza- 
tion of ig is such that the contraction of tg with four times the same antisymmetric 
matrix M yields 



tgMMMM := 24 tr M 4 - 6 (tr M 



2\2 



(42) 



Note a misprint in the definition of Z given in the appendix B of |26|: in Eq.(B.14) eio£io should 
be replaced by — eioeio. The sign which appears in the hep-th version of that paper is correct. We 
are grateful to P. Vanhove for a discussion on this point. 



15 



Eg = JpSiiEiiRRRR is shorthand for the so-called 'Euler-Lovelock density'; 
more explicitly, writing out all indices for once, we have 

E 8 = ±e ABCD ^ Ds e ABC E 1 ...E 8 RD 1 D 2 E ^---RD 7 D 8 E ^ 

= -SeI:::E & 8 Rd iD2 EiE2 ---Rd 7 d 8 E7Es m 

where 5| = +1, —1,0 is (when non-zero) the signature of the permutation be- 
tween D\ . . . D$ and Ex ... Eg. Observe that in Minkowskian signature the prod- 
uct £ii£n contains a minus sign with respect to the usual light-cone term entering 
string-theory amplitudes egeg = +<5f. 

There are many subtle issues related to the relative sign between the tgtgR^ = 
X and the \E% = —\Z contributions. In particular the usual lightcone calculation 
of 4-graviton amplitudes 11301 yields the same combination X — Z/8 for the IIA 
theory at tree level and at the one-string-loop level. However, according to Refs. 
Il39l l26l. the correct one-string-loop term (which lifts to the M-theory result) differs 
from the light cone calculation by an Euler-Lovelock term (and the Chern Simons 
term), and should read X - Z/8 + 2(Z/8) = X + Z/8 instead (this is precisely 
the opposite sign from the one obtained in the lightcone computation of [30]). To 
recover (or check) this result in the lightcone gauge one would have to calculate 
a 5-graviton amplitudes in IIA string theory at one loop, but we are not aware of 
such computations. This issue is important because, as we shall see below, our 
cosmological billiards are a priori sensitive to the Euler-Lovelock term, and hence 
to the flip of sign between the lightcone result X — Z/8, and the complete result 
X-Z/8 + 2(Z/8) = X + Z/8. 

In the following, we shall decompose the total contribution d39l into three sep- 
arate terms: the term L^ 1 = Jq = tstgR 4 + \E%, the Euler-Lovelock term 
£(4)2 = -Ie 8 = +2 ■ \Z contained in -2Z 2 , and the Chern-Simons term L cs 
contained in —2X2- The first term d40t (which coincides with the lightcone result 
X — Z/8) can be worked out more explicitly l39ll26ll4Tl . modulo terms that either 
vanish on-shell, or are equivalent to sub-leading terms (e.g. ~ E?F 2 ) 

= EJo = E (x - l -Z 
= VJLtj I — U L-AB ^CE ^DFGH 

+ AC abcd C a e c F Ce G b H C FG dh) (44) 

The second term C^ 2 has a much more complicated form when expressed in this 
way ll26*ll43l . and we will therefore leave it 'in the e-form' (177b . We shall discuss 
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the effects of this term in more detail in section 7, and of the Chern Simons term 
in section 8. 

The mixed terms containing (DF) factors have so far not been reduced to such 
a simple form. Let us just note that, in the (eleven-dimensional) light-cone gauge, 
the determination of all these terms can be reduced to an e trace over 50(9) V- 
matrices, and, more specifically, can be obtained by contracting any four factors 
from 

C a /3~,s := T^pT^Cijki , (DF) a/3l s := T^pTj™ DiFj k i m (45) 

with g ai — ai6 , as follows by inspection of the relevant vertex operators in 138*1 14T31 . 
Here the indices = 1, 9 label the nine transverse dimensions in eleven 

dimensions, and r iJ and T % ^ k are the standard 50(9) T-matrices with spinor in- 
dices a,@ = 1, . . . , 16. However, we would expect the caveats concerning the 
difference between light-cone and 'covariant' amplitudes, mentioned above for the 

4 contributions, also to apply to the mixed terms. 

5 Scaling behaviour of curvatures and field strengths 

Next we exhibit the leading asymptotic behaviour of the various fields and their 
'composites' corresponding to the higher order collections near T = 0, as func- 
tions of the logarithmic scale factors j3, making use of the decomposition (fl4l 
and other formulas derived in section 2. Exploiting the identification of the scale 
factors with the CSA degrees of freedom, we shall then try to interpret the expo- 
nential walls exp[— 2w(P)], as obtained from the higher order corrections, in terms 
of new exponential (Toda) walls exp[— 2a(j3)] in the E%q cr-model. To this aim, let 
us first determine the leading asymptotic behaviour of the curvature, starting with 
that of the spin connection uj and the anholonomicity coefficients Q,. Because the 
Iwasawa frame freezes in this limit, the 'leading order' amounts to neglecting the 
time derivatives of the intermediate frame, i.e. terms proportional to dt6 a . More 
precisely, the terms dt6 a are proportional to some exponential walls, starting with 
some simple root walls ~ exp(— 2aj(/3)). A consequence of this fact is that when 
discussing the subleading terms in the 8-th order correction (which are smaller than 

the dominant term by a factor ~ exp ( — 2 Wi<*t(/?)J with m > 0), one should, 
at some stage, take into account their 'mixing' with the terms proportional to dt9 a 
and coming from the leading curvature contribution. 

To proceed we now substitute the results derived in section 2 into the expres- 
sions for the Riemann tensor in order to see in more detail how its various compo- 
nents behave in the limit T — » 0. Retaining only the relevant terms containing time 
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derivatives, we get 



R-OaOb — —doUJabO ~ ^Oac^cbO + • • • — 

= 5 ab e 2a [d 2 t f3 a + d t ad t fi a - (d t f3 a f] + ■■■ (46) 



for the components with two O's, and 

R 0abc oc AT- V"^ 6 -^ = e'+P+P-F + permutations of a, b, c (47) 
for the component with one 0. For the purely spatial components we have 

Rabcd = Ra b l d + Rabid (48) 

where the first term on the r.h.s. contains the contributions with u a bo 

R abld '■= Uaco^bdo ~ Uado^bco = ^ c[a S b]d e 2a ' d t f3 a d t (3 b (49) 

while the second term R a b c d i nv °l ves ai l contributions from the purely spatial spin 
connection (1231 and is proportional to 

U abcd 01 e ° r 

OC ^+P b +W^-P d (5Q) 

together with all permutations of the indices a, b, c, d; note that some of these in- 
dices may coincide. The index e is independent of a, b, c, d; it is the summation 
index in the terms quadratic in the spin connection. For later reference, we note 
that the combinations appearing in the exponentials on the r.h.s. of (l50l differ from 
one another by linear combinations of 'symmetry roots' (f3 a — (3 h ), corresponding 
to the so-called 'symmetry walls' 0. 

As a check of our computation of the leading curvature components, let us 
work out the Einstein-Hilbert action to leading order; we get 



ER ~ e-^[Y,R { alb-^ R ^ 



a,b 



-2d 2 t a + Ys&n 2 - (d t a) 2 = -2d t a + £ G a bd t p a d t (5 h (51) 



a, b 



Therefore, modulo a total time derivative —2d 2 a, we recover the kinetic terms for 
the diagonal metric components in Q, with the usual (Lorentzian) 'superspace' 



18 



metric G a b from \2A\ . It is the inverse (= G ab ) of this metric which appears in 
both the Hamiltonians Q and © in the Introduction, and was written out in d2"oT ). 

It is easily seen that the leading components of the curvature are R a oaO an d 
R^bab ( w i tri a 7^ and no summation on repeated indices). These leading compo- 
nents are of order e 2cr and can be written as (remembering that, in our gauge, f3 a is 
approximately piecewise linear in the cordinate time t) 

RaOaO = e 2<T V a V a 

R abab = e2av a v b for a / & (52) 
where, for shortness of notation, we have defined 

v a :=d t p a , v a :=J2P b = dt<T-d t f3 a (53) 

b+a 

For the Kasner solution, the coordinate time can be taken to be t = — In T, and the 
'velocities' v a are then just the Kasner exponents. 

The remaining components of the curvature tensor are sub-leading in the sense 
that they are suppressed in comparison with the leading term (of order ~ e 2a ) by 
exponentially small factors. Before discussing in a unified way how one can define 
a relative ordering between several contributions, let us collect various results for 
the ratios between subleading curvature (and F-dependent) terms, and the leading 
curvature R ~ e 2cr . From d47b and d5Qb we see that 

Roabc I e 2a ~ e -°+P a +P b -P c + permutations (54) 

Notice that the combination appearing in the exponential on the r.h.s. is just (mi- 
nus) the gravitational wall form d35b . Similar estimates hold for R^ b \ d \ a represen- 
tative example is 

/ ^ ~ e-^+P+P+P-P 1 (55) 

and all other such quotients differing only by symmetry walls. 

The asymptotic behaviour of the 4-form and of its derivatives can be analyzed 
in a similar manner. For the electric field strength we have 

Foabc = e^ a ~^ c F tabc (56) 

where 

Ftabc '■= QmaOnbOpcFt mnP (57) 

and the overall factor e a comes from the inverse lapse N~ l . [We recall that the in- 
dices 0, a, b, c denote orfhonormal-frame indices.] As shown in [5 1, it is precisely 
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the barred field strength F ta bc (in the gauge N = yfg; or, in any gauge, the con- 
jugate momenta N^fg F tmnp ) which possesses a limit as T —> 0. Comparing the 
electric field strength with the square root of the leading curvature R 1 / 2 ~ e a , we 
obtain [no summation on (a, b, c)] 

F 0abc I e° ~ e -V+?+fn (58) 

where the combination (P a + (3 b + (3 C ) appearing in the exponential on the r.h.s. is 
related to an 'electric wall', see 0. In particular, for (a,b,c) = (1,2,3) we just 
re-obtain the simple root ao in d28t . which explains the term 'electric root'. We 
thus conclude that the electric field strength is 'smaller' than the square root of the 
dominant curvature by one electric wall. 

For the magnetic field strength, we have similarly 

Fabcd = <P + * + P+*F abei (59) 

with 

F a bcd '■= Ga m 9 b n 9c P &d q Fmnpq (60) 

Again the barred field strength is the one that remains finite on the initial hyper- 
surface. Note the different positions of the world indices in these two definitions. 
[The correct index position for the quantities that freeze is easy to remember as the 
quantities that freeze are the Q's and the P's: A mnp and their spatial derivatives, 
and their conjugate momenta Ti mnp .] Comparing the magnetic field strength to the 
square root of the leading curvature, we obtain 

Fabcd I e a ~ e -S^a,b,c,<i^ (61) 

The combination appearing in the exponential on the r.h.s. is now a magnetic wall, 
which is also equal to the sum of two (non-overlapping) electric walls. Conse- 
quently, the magnetic field strength is 'smaller' than the electric field strength by 
one electric wall, and 'smaller' than the square root of the leading curvature term 
by two electric walls. This shows that all the higher order terms which are algebraic 
in R and F are sub-dominant compared to the leading R A terms. For instance, F 8 
is 'smaller' than the leading R 4 term by eight electric walls. 

Finally, the covariant derivative of the field strength is dominated by the time 
derivative of an electric component: 

D F 0abc = e 2a - pa -P b -P c D t F tabc + ... (62) 

Hence 

D F 0abc I e 2 ° ~ e-C^+zn (63) 
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In other words, DF is smaller than R by one electric wall. 

We can conveniently summarize the relative ordering of the building blocks en- 
tering the 8th order correction by introducing the symbol ■< to denote the ordering 
of exponential factors involving simple wall forms. As will be further discussed 
at the beginning of Section 6, we can treat the negative exponentials of the sim- 
ple roots as small parameters: £j ~ exp(— oti{(3)) < 1. This defines an ordering 
such that m > •n! i implies exp(— ^ mai(f3)) -< exp(—J2 n i a i(P))- Now, by 
looking more closely at the various exponential factors appearing on the r.h.s's of 
the comparison equations <54b . d55b . d58t . doTT) . d63l above, one sees that they can 
be rewritten as negative exponentials of sums of simple roots. For instance, the 
r.h.s. of (154} is e -(°-P a -P b +P c ) ^ where the combination a - (i a - (3 b + (3 C is the 
gravitational wall form, see d35l above, which can be written as a sum of simple 
roots, see Eq. J36t , Actually, it is convenient, for a quick perusal of the relative 
ordering of various contributions, to neglect the ordering with respect to symmetry 
roots ai,...,ag, see d29t in the next Section, and to focus on the ordering defined 
by the number no of leading electric roots ao(P) = j3 l + 1 + f3 3 . In other words, 
we shall simply write exp(— Yl n i a i(fl)) di exp(— noa>o(f3)), and use the order- 
ing exp(— noao(P)) d exp(— n' ao((3)) when n' > uq. For instance, in view of 
the fact that the gravitational wall form which appears in (I54t is at level t = 3 
(see d36l ). we can write e~( CT-/3<1 ~^ + P C > ^ e~ 3a °. Finally, with this notation, the 
above results simplify to 

Roabc I e 2a 1 e- 3a ° (64) 

R [ ah I ^ < e~ 6a ° (65) 
Fl bc /e 2 °<e- 2 ^ , Fl bcd /e^<e~^ (66) 

(in comparing with the curvature which is second order in derivatives, we must 
square the 4-form field strengths because they have only one derivative) 

DF I e 2a < e~ a ° (67) 

These estimates show that the dominant terms among the various 8th-order 
corrections R 4 , R 2 (DF) 2 , R(DF) 3 , (DF) 4 , ■ ■ ■ , F 8 are the curvature terms R 4 . 
Among the R 4 terms, the leading contributions will come either from the time- 
time components Ro a Oa, or from the purely spatial ones ~ ^ilL,- Then, the next 
to leading contributions will come either from R?(DF) 2 or from R?F 2 , and will 
be smaller than the leading one by a factor ~ e~ 2a °. 
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6 Leading behavior at higher orders: curvature terms 



Before proceeding with the analysis of the higher order derivative contributions in 
the Lie algebraic setting, let us clarify in which regime our analysis of their ef- 
fects, and in particular of the 8-th order derivative corrections, will be meaningful. 
This regime is one of an 'intermediate asymptotics' of the type Tp <C T <C To- 
Here, Tp is the (eleven-dimensional) Planck time, and To is some initial (or final) 
time away from the big crunch (or big bang) singularity located at T = 0. We 
assume that the Cauchy data at To are such that spatial derivatives of the metric 
and the 3-form are mildly smaller than their time derivatives, say d a g < 0.2drg. 
We also assume that the time derivatives correspond to curvatures much smaller 
than the eleven-dimensional Planck curvature T p . Then, for a while, the BKL- 
type analysis becomes better and better for times T <C To, in the sense that the 
potential walls become sharper and sharper, and that the Weyl chamber of E\q 
appears as a dominant billiard description of the dynamics of the diagonal com- 
ponents of the metric at any given spatial point. Such a billiard description of 
the Einstein-three-form system remains accurate until the higher-order corrections 
to the action, ~ R + Tpi? 4 + • • • become important. This happens for curva- 
tures R ~ Tp 2 , i.e. for cosmological times T ~ Tp. As there are no small 
coupling parameters in M-theory, when the leading correction ~ R 4 becomes im- 
portant, all higher-order ones ~ Tp N ~^ R N become also important as T ~ Tp. 
However, in the intermediate regime Tp <C T <C To, it is meaningful to con- 
sider only the terms ~ R A and to treat them as a small correction to the usual 
supergravity billiard defined by So = R + F 2 + AFF. In other words, in such 
a regime we can meaningfully combine two separate expansion schemes: (i) a 
height expansion in which the negative exponentials of the simple roots are treated 
as small parameters, say £j ~ exp(— cti{0)), and (ii) a small curvature expan- 
sion in which e c ~ TpR ~ (Tp/T) 2 is treated as a separate small parameter. 
Note that the basic rule behind the first expansion is to order the exponential 
terms by their height. Indeed, a generic wall exp(— 2a (/?)), with a positive root 
a = niCti is seen as the product of small numbers, nfe) 2 ™ 1 ' which gets smaller 
and smaller as the rij increase. Another, and more heuristic, way to think about 
the height expansion would be to replace the ten individual (3 a, s, or better their 
simple-root combinations aj(/3), by some 'average scale factor' (3, or some 'aver- 
age simple root' a s = a-j. This is a legitimate approximation because the oti(f3) 
increase toward +oo and remain positive during their chaotic motion within the 
Weyl chamber a« {0) > 0. In terms of such averaged quantities one can estimate 
the effect of a generic wall exp(— 2 ^ rajO!j(/?)) as being on average equivalent to 
exp(-2(X>i)"s) = exp(-2ht(a)a s ). 

As we have seen in Section 4, the dominant terms, in the sense of the height ex- 
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pansion just explained, among the 8th order derivative corrections R , R 2 (DF) 2 , 
R(DF) 3 , (DF) A , ■ ■ ■ , F 8 are the curvature terms i? 4 . Therefore, these are the 
terms that will define the leading corrections to the dynamics in the intermediate 
asymptotics T P <C T <C T . 

Let us examine more carefully the contribution of the terms R A in the Hamilto- 
nian constraint TL. In any gauge an additional term J in the (invariant) Lagrangian 
density contributes the term —EJ in the Hamiltonian constraint. We recall from 
( flTt that the eleven-dimensional volume density behaves as E = N \Jg = e~ 2cr . 
As we have seen above, the leading frame components of the curvature blow up 
proportionally to e 2a '. Hence, a term J ~ R A contributes —EJ ~ e~ 2cr (e 2cr ) 4 ~ 
e 6<J . [Certain special quartic curvature invariants R 4 might only contribute terms 
<C (RoaOa) 4 , but we shall see below that both the tgtgR 4 and £§ do contain terms 
~ (e 2cr ) 4 .] For more generality, let us consider a correction containing an arbitrary 
power of the curvature, say J ~ R (the precise kinematical structure of R does 
not matter at this point). This contributes to the Hamiltonian constraint a term 

ER N oc e 2 ^- 1 ^ (68) 

The important point for our analysis is that, in the Hamiltonian constraint such 
a term has the form of an exponential wall ~ e~ 2w<y ^\ where w {(3) is a linear form 
in the /3's. The associated 'wall form' wn = w[R N ] is easily read off, viz. 

w N ((3) = -(N-l)a <{=> w N = -(N- 1)- (1,1, 1,1, 1,1, 1,1, 1,1) (69) 

Let us now see for which values of N the wall form d69l happens to be a root of 
Eiq. A simple necessary criterion for checking whether a wall form w might be a 
root of Eiq consists in computing the squared length of w with the scalar product 
of the CSA of E 10 (|26l>. Indeed, all forms w on the CSA which belong to the 
.Eio root lattice 7 , and which are such that w 2 = —2j, for some integer j > 0, are 
(imaginary) roots of E\q and vice versa. For the wall form d69t . one easily finds 

*4 = -y(iv-i) 2 (70) 

This is negative, and is an even integer only when N = 3k + 1, with some positive 
integer k. This includes the case of main interest here, namely N = 4, i.e. k = 1. 

It is easily checked that when the necessary condition N = 3k + 1 is met, the 
corresponding w^k+i = —3ka does belong to the root lattice, and is therefore an 

7 We recall that the root lattice, i.e. the set of integer combinations of simple roots, contains, but 
is strictly larger than, the set of roots. 
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(imaginary) root. For instance, the quartic correction R 4 corresponds to the Eiq 
root 



which is of squared length a 2 = —10 and height ht(a) = —115; its level is 
i = — 10. We see that this root lies rather deep inside the lightcone in root space. 
More generally, for N = 3k + 1 (i.e. for corrections R , R , R 10 , ■ ■ •), one gets 
wall forms corresponding to E\q roots of squared length — 10k 2 , height ht(a) = 
-115/c, and level t = -Wk for k = 1, 2, 3, . . .. 

Looking back at the potential wall do^t associated to the correction R N we 
notice something special. Barring special cancellations, all terms in the correc- 
tion R N lead to a unique dominant potential wall e 2 ( N ~ l ) u . This is quite dif- 
ferent from what happens for the usual walls that have appeared in the analysis 
of [8 5|. Indeed, the walls associated to the electric energy of the 3-form, or 
its magnetic energy, or to the 'gravitational energy' all came in multiplets of the 
group of permutations of spatial indices. For instance, the electric walls were of 
type exp(-2(/3 a + [3 b + (3 C )), where a, b, c are three different spatial indices. The 
lowest-height electric root a>o = (3 1 + (3 2 + (3 3 is accompanied by other electric 
roots at the same level I = 1, which are obtained by acting on oq by the generators 
of SX(10) (with roots (3 2 - (3 1 , ■ ■ •), which have the effect of permuting the SX(10) 
indices among a = 1, 2, • • • , 10. This fact was associated to the fact that the basic 
electric root o.q = j3 l + (3 2 + /3 3 was simply the lowest weight associated to a 
generator of Eiq which was given by an antisymmetric 3-tensor 5X(10) represen- 
tation, say E abc . In other words, the walls which come in permutation multiplets 
are associated to non-trivial tensor representations of SL( 10). 

The 'isotropy' of the wall form d69l . i.e. the fact that it does not depend on 
the choice of a special permutation of SX(10) indices, suggests that wn, Eq. d69l . 
should be associated to a root of Eiq whose associated generator is a SX(10) sin- 
glet. We do find in the tabulated results of l24ll that there are indeed the requisite 
singlets matching the R 4 terms at level I = 10, and the predicted R 7 terms at level 
i = 20. Moreover, these are the only singlets among the 4 400 752 653 represen- 
tations identified up to level I < 28 in [24 1. In view of the 'contravariant' char- 
acter of the 'raising generators' of E w , namely, E^ = E abc , E^ = E ai - a6 ~ 
[E( 3 \ E^j, etc. it is easy to see that SX(10) singlets can only occur at levels 
£ = 10/c (so that one can use 3k covariant e tensors of SX(10) to soak up the 3£ 
contravariant indices of E^) 8 . It is, however, remarkable that: (i) singlets do 
occur at these levels with non zero outer multiplicities, and (ii) there exists such 

8 We thank Axel Kleinschmidt for pointing out this algebraic fact. 




(71) 
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a nice compatibility between the algebraic structure of E\q (built by 'commuting' 
three-form level- 1 generators E abc ), and the wall structure entailed by curvature 
corrections. See the final Section for further discussion of this point. 

Moreover, among the highest weights of the representations at levels i = 10 
and t = 20 the associated root (TTTb and its integer multiples are the ones of largest 
root height. Therefore, for any N, among the roots associated with the 2N-th 
order derivative corrections, ( 1771) is always the root that lies deepest inside the 
past lightcone. For the moment, no data are available beyond I > 28, but we 
confidently predict the requisite 5L(10) singlets to appear at levels lOfc with roots 
3ka and non-vanishing outer multiplicities. 

We can summarize our findings in a little lemma. 

Lemma: Corrections proportional to the N-th power of the Riemann (or Weyl) 
tensor are compatible with E\q if and only if N = 3k + 1, i.e. for R 4 , R 7 , R 10 , . . .. 
Moreover, such corrections must be associated with sl\o-singlet generators of E\q 
at levels i = — 10k. 

For the special value k = 0, the 'leading term' is just the kinetic term d5Tl) . and 
there is no 'leading wall'. For k = 1, 2, i.e. N = 4, 7, the outer multiplicities of 
the singlets are known: fi = 3 for £ = 10, and fi = 913 for i = 20 |24|. We have 
no physical interpretation for these values in terms of the R 4 and R 7 corrections to 
the action. It might have been tempting to associate the outer multiplicity \x = 3 to 
the three different contributions (tgig-R 4 , E$ and the Chern-Simons term) entering 
d39b . However, we shall see that the Chern-Simons one is subdominant. Note 
also that for E\\ the root corresponding to (TTTb is associated with the SX(ll) 
representation (0000000003), which is no longer a singlet. 

Below, we will show that this result extends to corrections of mixed type 
R N ~ n (DF) n . Remarkably, the very same constraint N = 3k + 1 on the allowed 
powers of curvature, which we found here assuming compatibility with the Eiq 
root lattice, was arrived at by a very different argument in [ 47 1 (to wit, an expan- 
sion in integer powers of the inverse membrane tension). Note also that, as far 
as we know, no detailed information is yet available about the R N corrections for 
N > 4. 

7 Positivity properties of leading curvature terms 

Let us comment on two interesting aspects of the Ei root dTTl associated with 
the leading R 4 corrections. First, as announced in the Introduction the root (17 It . 
besides being imaginary (or more precisely, time-like, a 2 < 0), is also negative, 
a G A_. Physically, this corresponds to the fact that the R 4 corrections grow 
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faster than the leading Einstein terms near the singularity. As we already explained, 
this fact limits the validity of our analysis to some intermediate asymptotic Tp <C 
T <C Tq. Within such an intermediate asymptotic expansion, we see no a priori 
conflict between the occurrence of negative E\q roots in the M-theory action, and 
the form © of the coset action, involving only positive roots of _E 10 . Indeed, as 
already pointed out above, the positive root form © is the result of using a special 
nilpotent (i.e. purely upper triangular) parametrization of the coset E\q/ K(Eiq). 
The use of a more complicated parametrization, involving a mixture of upper and 
lower triangular parametrization, would entail the presence of negative roots in the 
coset action. 

However, no matter how one parametrizes the solution, the trajectory that ex- 
tremizes the coset action is simply a null geodesic on E\q/ K{E\q), and therefore 
geometrically the same solution. In the case of finite-dimensional Lie algebras one 
can prove that a generic geodesic curve is conjugate to one in the CSA. However, 
no proof of this result is known for hyperbolic Kac-Moody algebras, nor is it known 
whether such a statement holds at all, see [48 1. If that were true for geodesies on 
E\q/ 'K(Eiq) one would immediately conclude that a generic null geodesic is con- 
jugate to a null geodesic in the CSA, i.e. to the curve j3 a (t) = v a t + f5 a °^ (with 
v 2 = 0) which runs straight from 'past null infinity' towards 'future null infin- 
ity' in the CSA. Fortunately, we can qualitatively control the global behaviour of 
null geodesies without having to assume the unproven equivalence with a CSA 
geodesic. 9 Indeed, we can remark that, in the Borel parametrization of Eq. ©, 
the potential terms oc II 2 are all non-negative, so that they define a non-negative 
'squared-mass term', — tt 2 = M 2 > for the motion projected in /3-space. As 
f3 a oc ir a this implies that the f3 a worldline is everywhere timelike or null. Re- 
membering that, in terms of the logarithmic scale factors /3 a 's, 'future infinity' is 
a = Pa — ► +oo, and corresponds to a 'big crunch' singularity where the vol- 
ume Jg = e~ a — ► 0, we conclude that a generic null geodesic, which is initially 
future-directed, will never be able to turn back to become past-directed. In other 
words, this proves that a generic E\q null geodesic cannot correspond to a 'bounc- 
ing' cosmological solution, where the volume of space reaches a minimum before 
it re-expands. In view of this argument, and if there is any truth to the E\q con- 
jecture, we expect that the sign of the higher-order corrections to the supergravity 
action will be such as to allow no bouncing solutions. 

We have studied the sign and magnitude of the various leading quartic correc- 



9 See, in this respect, the Appendix below where we emphasize that null geodesies on 
Ew/K(Ew) are 'strongly chaotic' in the sense of being extremely sensitive to small changes in 
their initial conditions. This result suggests that geodesies which do not lie within the CSA might be 
somewhat different from the simpler (Kasner) CSA ones. 
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tions in (l39l to see if they might be compatible with our expectations. 10 Actually, 
the analysis of the sign of the quartic-in-i? contribution to the gravity action is 
rather involved, because the precise kinematical structure of the R 4 now does mat- 
ter. The first important fact is that the sign of the coefficient of the loop correction 
is positive 1 1 . Then there arises the issue of the actual sign of the two leading 
contributions to the action, d40b and the Eg term in d4lT) (which, as we shall see 
below, is dominant compared to the Chern-Simons one) when they are evaluated 
'on-shell', i.e. along a solution of the leading Einstein-3-form dynamics. If we can 
prove that the sum of the on-shell actions d40b and the Eg term in d4lT) is numeri- 
cally positive, this will prove that the corresponding contribution to the correction 
of the gravity Hamiltonian is 12 negative. In that case, the wall associated to the R 4 
corrections (characterized by a negative, time-like root, with negative coefficient) 
will represent no longer a 'mountain', which would obstruct the motion of to- 
wards infinity and thereby induce a bounce, but instead a 'canyon' or a 'crevice' 
precipitating the collapse of space toward zero volume (i.e. a — > +oo). 

Let us therefore have a closer look at the various R 4 contributions to the action. 
We consider first the contributions coming from (140b . which explicitly read as d44b . 
and will later consider the E$ one contained in (EiTli . The first observation is that 
in leading order (in the sense explained above), only terms with i? a 0a0 an d R^bab 
need to be considered in d44t . Secondly, it is straightforward to verify that all terms 
proportional to either (Co a oa) 4 or (C a b a b) A (no summation on indices!) cancel 
between the two contributions in (144b . Substitution of d52b into d44b and some 
further calculation then reveals for the leading part of the density d40t the explicit 
expression 

EJ = 768e 6CT T/ 2 4 , (72) 

where 

V 4 = ^ v l v b( v aVa + V b V b ~ V a V b ) 2 + ^^2 V a v b v c( v a + V b + V c ) 2 > (73) 
a, b a,b,c 

and the sum ^ over a, b is restricted to a ^ b, while the sum J2" ranges over the 
three indices a, b, c with the restriction that a ^ b,b ^ c, and c ^ a. Note that each 

'"Strictly speaking, the sign of the leading correction is not logically related to the bounce/non- 
bounce behaviour. Indeed, even if the leading correction indicates a bounce behaviour, the next order 
corrections may prevent the completion of this bounce. 

"Due to a remnant Euclidean-type convention, one should reverse the overall sign of the actions 
in (39). 

12 We use here the usual first order result Hi — — Li for the correction to the Hamiltonian cor- 
responding to a correction to the Lagrangian. This result applies to the case where one introduces 
no new degrees of freedom, but only modifies the action for given degrees of freedom, possibly by 
higher derivative terms. 
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ordered pair a < b contributes two equivalent terms to the sum ^ , while each 
ordered triplet a < b < c contributes six equivalent terms to the sum ]P . The 
explicit expression d73t as a sum of squares proves that Jo > 0. If that were the 
only contribution to the action, this argument would suffice to prove the negativity 
of the corresponding term in the Hamiltonian. However, we still need to consider 
the term proportional to the Euler-Lovelock density in (14 It . Using the explicit 
expression of the Euler-Lovelock densities in Kasner spacetimes [49], 13 we get the 
following explicit, leading expression for the relevant term in PTt 



£( 4 ) 2 = --EEs = +^EZ = -46080 [d t (e 6a V 7 ) + e 6a V 8 ] (74) 

where 

V 7 := ^2 v ax---v a7 (75) 

ai<---<ay 

and 

Vs ■= ^ v ai ■ ■ ■ v as (76) 

ai<-"<08 

Note that, contrary to what happened for the sums ^ and above, the sums 
appearing in Vj and Vs contain only one contribution per ordered multiplet of in- 
dices. In the following, we shall discard the total time derivative 14 appearing as 
first term on the right-hand side of (1741 . and evaluate 

£ (4)2' = _4 60 80e 6CT V8 (77) 

In order to see whether the term ( 1771 is likely to modify the positivity of (1731 
we need to compare (46080/768) Vs = 60 Vs to the term V^ 4 , i.e. to the expression 
defined by the restricted sums in dT^l . The latter quantity, V^ 4 , was shown 

above to be a sum of squares, and therefore to be always > 0. We need to study 
the sign and magnitude of the additional contribution —60 Vs. 

We can obtain rather simply an approximate knowledge of the magnitude of 
Vs in the following way. 15 We consider the 10-th degree polynomial P(x) = 
Q a (x — v a ). Let us consider on-shell values of the v a 's, i.e. such that v 2 = 
— 2J2 a <b v a v b = 0- Let us also normalize the values of the v's by imposing 

13 After discussion with the author, it appears that the overall coefficient 2 2p ~ 1 in Eqs. (22a) and 
(22b) of 1 49] should be replaced by the, alas, larger coefficient (2p)!/(2p — 1). 

14 Note, however, in passing that if one were to evaluate £S^ 2 on-shell before discarding this term, 
it would yield a non-zero contribution (but a contribution which is equivalent to a reparametrization 
of the basic degrees of freedom). 

15 We thank Ofer Gabber for suggesting this method for bracketing the value of Vs. 
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Y^ a v a = 1- [So that we have both Yla v a = 1 = J2a v a}- m tne expansion of 
the polynomial P(x) = x 10 — cgx 9 + cgx 8 + • • • + C2X 2 — c\x + cq, we then 
know c 9 = J2a v a = 1. c 8 = Sa<6 ^6 = and we have c 2 = V 8 , where y 8 
is the quantity we wish to bracket. Because all the roots of P(x) are real, it is 
well known that all the roots of its successive derivatives P'(x),P"(x), ■ • • will 
also be real. By considering the expansion in powers of x of the second derivative 
P"(x)/90 = x 8 - 0.8x 7 + Ox 6 + • • • + Vs/45, and relating it to its real roots, 
say Xi (for i = 1, • • • , 8), we get the following expression for Vg: Vg = 45 f\i %i 
where the Xj's are constrained by ^ Xj = 0.8 and J2i x i = (0.8) 2 . It is then easy 
to get both upper and lower limits to Vg by considering the various possible cases 
where p values among the Xj's are negative. [It is easily seen that p must be p > 1.] 
Using then the convexity inequality a\a^ ■ ■ ■ c? n < [(Xw a i)/ n ] n -> w i tn n = 8 — p 
and/or n = p, and minimizing over the average value of the negative Xj's, we get 
the double-sided inequality 

- 45(0.8) 8 /8 4 - 5 < V 8 /(a v ) 8 < +45(0.8) 8 /8 5 (78) 

where we denote a v = ^ a v a- Numerically, this yields for the quantity we are 
interested in: 

- 0.0391 < 60V 8 /(a„) 8 < +0.0138 (79) 

These rather small allowed values for V% are due to the fact that V% is a sum of 
products of 8 different u's. By contrast, we expect the quantity y 2 4 / (ct„) 8 in d73b to 
typically take significantly larger values because it contains products of only two 
or three different v's. This suggestive (if not quite conclusive) argument makes 
it plausible that the difference y 2 4 — QOVs will be always positive. Actually, the 
place where this difference is most likely to become negative is in the vicinity of 
the points where (V2) 4 vanishes. This happens (normalizing again by ^ a v a = 1) 
when one specific v is very close to 1 , and the other ones are near zero. [Note that 
the Kasner solution where (v a ) = (1, 0, 0, • • • , 0) is flat spacetime in disguise.] 
However, it is easy to see that, near such a point, V 2 4 = 0(v 4 ) while Vg = 0(v 7 ). 
Therefore, near the points where K 4 is the smallest, the difference V^ 4 — 60 Vs 
remains positive. We performed some partial numerical investigations of the sign 
of the difference V^ 4 — 60 Vg and always found it to be positive. 

We take all this as evidence that the total 8th order action d39l is equivalent to 
a negative contribution in the Hamiltonian, i.e. to an 'inverted wall', or a 'crevice'. 
As we said above, this is compatible with the E\$ expectation that there should be 
no bounce. 
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8 Subleading terms: some examples 

We now wish to rephrase the results derived in section 4 in a Lie algebraic setting, 
by giving more details concerning the subleading 8th-order derivative corrections, 
and by illustrating their asymptotic behavior with some examples. More specifi- 
cally, we will demonstrate explicitly how all these terms can be associated either 
with (in general imaginary) roots of E\q, or, more generally, with elements of the 
E\q root lattice. Comparing the relevant highest roots with the corresponding high- 
est weights in the tables of |24| we can then read off to which s[\q multiplets they 
belong when they do correspond to actual roots. 

In section 4, we estimated the various tensor structures in comparison with the 
leading curvature terms R ~ e 2<T . We have already remarked that the subleading 
terms are smaller than the latter terms by exponential factors involving electric 
and symmetry wall forms, which are linear combinations of the simple wall forms 
appearing in (l29l . Above, we discussed a coarse ordering of the effects of these 
wall forms by focussing on the number of basic electric roots «o entering them. We 
shall now have a closer look at the various wall forms associated with the different 
types of higher-order corrections, and their relation to the root lattice of E\q. 

From the results in Section 4 above, and from the basic fact that in the final 
expression all indices must be contracted in pairs, and thus appear an even num- 
ber of times, we can see that the subleading terms in the Hamiltonian constraint 
are all obtained by multiplying the leading Hamiltonian contribution e 6fJ (corre- 
sponding to the leading root ai ea ding = — 3d) by an exponential factor of the form 
exp(— 2 n%oii), with nj being some non-negative integers. Because of the pres- 
ence of a factor two in the exponent (linked to the fact that all indices are finally 
contracted), we see that all the wall forms associated to the corrections can be 
written as exp(— 2u; su bi ea di ng ), with subleading wall forms of the type 



We see from (l8Ql that the set of wall forms associated with higher-order corrections 
is a subset of the root lattice of E\q. To better understand this subset, it is conve- 
nient to introduce at this point a geometric picture of this subset. This is done in 
Figure 1. 

The reasoning above shows only that the wall forms are of the form (Op for 
some non-negative integers n%. The tools we used above are too coarse to de- 
termine whether all values of the integers ni are realized, or whether there are 
restrictions of any kind on the values of the n^'s. In view of the complicated struc- 
ture of the higher-order corrections, which as we pointed out, is not known in an 
algebraically completely reduced form, especially for the crucial next-to-leading 




(80) 
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Figure 1: Schematic representation of the roots of Eio in the Lorentzian space of the 
Cartan degrees of freedom j3. Positive roots are in the upper part of the figure, while 
negative ones are in the lower part. The real roots are on the outer (timelike) hyperboloid 
a 2 = +2. The imaginary roots are either null {a 2 = 0) or timelike {a 2 = — 2j). The basic 
root associated with the leading R A correction to the supergravity action is the negative, 
timelike point — 3er indicated. From it stems the solid half-hyperboloid of roots associated 
with subleading contributions contained in R A as well as in the other 8th order derivative 
terms: R?{DF) 2 , R(DF) 3 , (DF) 4 , ■ ■ ■ ,F S , as well as the Chern-Simons contribution 
AR A . Similar half-hyperboloids stem from the roots — 3ka predicted to exist and to be 
associated to higher-order contributions ~ R 7 , R 10 , ■ ■ ■ 
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terms R 2 (DF) 2 and R 3 F 2 ), we have made no attempt at an exhaustive analysis 
as to whether the integers n.j are indeed subject to restrictions. Instead, we have 
proceeded 'experimentally' by studying in detail some specific, but typical, terms 
that we could check to be indeed present among the available expressions for var- 
ious higher-order corrections. In all the cases we have checked, we did find a 
rather remarkable pattern: the vector in root space joining the 'leading ' i? 4 root 
^leading = — 3cr to any subleading higher-order wall form u> su bieading, i-e. the 
quantity 



in all cases. Below, we will refer to ( as the 'relative vector'. This leads us to 
formulate the following 

Conjecture: The relation ( 182 I I holds for all kinematically allowed combinations, 
and thereby encodes information about the algebraic structure of the 8th order 
corrections. 

We shall illustrate below with some examples the interrelation between this 
conjecture and the properties of the E\q root lattice on the one hand, and certain 
'kinematical cancellations', as they follow from the kinematical structure of the 
known 8th order correction terms, on the other hand. However, we will leave a 
systematic investigation of this aspect to future work. 

Geometrically, the restriction ( l82t (together with n« > 0) means that the rela- 
tive position vector (, Eq. d8Tb . belongs to the set of positive roots, A + . In other 
words, if we assume that this restriction is indeed true for all terms, we can geo- 
metrically describe the set of wall forms generated by higher-order corrections as 
a solid half-hyperboloid, which is congruent to A + , and with the leading R 4 root 
^leading = — 3(7 as its basis, see Figure 1. From this picture there seems to be 
no evident upper limit on the height ^ n-i of the relative vector (UTTl . We initially 
thought that the wall forms ^subleading might be constrained to lie inside the root 
diagram stricto sensu (instead of lying anywhere in the root lattice), which would 
have implied the further restriction (u> su bieading) 2 < 2. However, as we shall see 
below, we found wall forms lying high enough in the solid half-hyperboloid d82t 
to extend beyond the outer hyperboloid w 2 = 2 corresponding to real roots. More 
precisely, by adding simple roots of Aq one generates slio weight diagrams start- 
ing at a certain highest weight inside this solid half-hyperboloid. Many of the slio 
weight diagrams are such that only a subset of the weights are actually roots of 




(81) 



was 'experimentally' found to satisfy the restriction 




(82) 
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Eio, while the weight diagram itself extends beyond the hyperboloid w 2 = 2. 
We emphasize that there is no a priori inconsistency or incompatibility in this fea- 
ture: indeed, starting from a Hamiltonian containing only roots as wall forms, any 
change of parametrization can generate arbitrary combinations of the roots, and 
therefore arbitrary elements of the root lattice. 

Let us now substantiate our conjecture ( 182 1 by giving examples of wall forms 
associated by various higher-order contributions. We first consider the Chern Si- 
mons term in (PUT) , which is proportional to 

ABCDl.-Ds a n E 1 E 2 p E 7 E 8 .(S) 

= e u A 3 [24 tr R 4 - 6 (tr R 2 ) 2 ] (83) 

Here the tensor tg is, by virtue of d4*2"l) . equivalent to the combination of traces ap- 
pearing in d4lT) . Though the above contribution is also quartic in the curvature, and 
contains the spatial components A a t, c of the three-form which simply freeze near 
the singularity, it is actually subleading w.r.t. the leading curvature contribution 
~ ER 4 ~ e 6cr . The reason for this is the interplay of the e tensor, which constrains 
the first pairs of indices [D1D2] • • • [D-fDg] on the curvature tensors to be all differ- 
ent, with the effect of the tg tensor, which contains Kronecker deltas and thereby 
obliges the remaining indices on the curvature tensors to coincide pairwise. As a 
consequence, it is not possible to have only leading curvature components of the 
type RoaOa or R a bab- A typical contribution is 16 

(84) 

where we used <47b . After reshuffling some indices, we obtain the wall form 

w = -(3, 2, 2, 2, 1,1, 1,1, 1,1) w = -[5; 2,5,8,6,5,4,3,2,1] (85) 

This is a root because w 2 = 2, with I = —5. The corresponding relative vector 

C = w- ai eadin g = w + 3a = (0, 1, 1, 1, 2, 2, 2, 2, 2, 2) (86) 

satisfies (w + 3a) 2 = +2, in agreement with d82b . We thus see that this Chern 
Simons wall form is so much higher than the basic R 4 root —3a that the associated 
wall form lies even beyond the interior of the past CSA lightcone and belongs to 
the timelike hyperboloid a 2 = +2. 

Next we consider the subleading curvature terms. Let us first recall that the 
kinematic structure of the C 4 terms in d44l) eliminates certain combinations; for 



6 In the remainder, we will no longer distinguish between the Riemann and the Weyl tensor. 
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instance, from d73l it is obvious that there are no contributions of the type C^ bab (no 
summation on indices a, b). As is evident from d50t . the curvature terms coming 

(2) 

from R abcd and the above estimates are strongly suppressed. Furthermore, not all 
of them produce actual roots, as we shall see presently. 

Obviously, from (BOl . there are many such terms, which we can group into 
certain 'permutation multiplets'. A first set of terms corresponds to terms of type 

R d(2) n(2) d(2) „(2) 2 <r . /3 a +/3 6 +/3 c -/3 d . /3 a +/3 f '+/3' i + / 3/-2/3' 

We can write this result schematically as 

4 



^ e -2a . e 2(/3"+/3''+/3 c +/3 e +/3-f+ / 39+/3' l - / 3 l - / 3J-/3 ft ) (gg-j 



Taking the indices to be all different, and choosing k) = (8, 9, 10), this gives 
w = -(0,0,0,0,0,0,0,2,2,2) w = -[2; 2,4,6,6,6,6,6,4,2] (89) 

which is clearly not a root as w 2 = +8. Nevertheless, the relative vector 

£ = w + 3a = (3, 3, 3, 3, 3, 3, 3, 1, 1, 1) (90) 

obeys ( 2 = 2 and therefore lies on the boundary of the half-hyperboloid in agree- 
ment with our conjecture. A second set of terms derives from 



E 



-2a . +2(p a +p b +p c +p d +p e +pf~pa-P h ) 



(91) 



and differs from the one above by one symmetry wall. Again taking all indices 
different and choosing them appropriately, we get 

w = -(0,0,0,0,0,0,1,1,2,2) w = -[2; 2,4,5,6,6,6,6,4,2] (92) 

which is again not a root as w 2 = +6. The relative vector 

£ = w + 3o- = (3, 3, 3, 3, 3, 3, 2, 2, 1, 1) (93) 

is lightlike, Q 2 = 0, and lies within the half-hyperboloid. A third set of terms 
comes from 



E 



^ e -2a . e +2(/3 a +/3 b +/3 c +/3 e +/3 9 -/3' 1 ) 



(94) 



and gives rise to, choosing values for the indices conveniently, 

w = -(0,0,0,0,0,1,1,1,1,2) w = -[2; 2,4,6,6,6,5,4,3,2] (95) 
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with w 2 = +4, and 

C= (3,3,3,3,3,2,2,2,2,1) (96) 
with ( 2 = —2. Finally, there is a fourth set with 

w = -(0,0, 0,0, 1,1, 1,1, 1,1) w = - [2; 2,4,6,6,5,4,3,2,1] (97) 

which now is a root because w 2 = 2, and 

C = (3,3,3,3,2,2,2,2,2,2) Q 2 = -A (98) 

again confirming the conjecture. The maximal height for any of these wall forms 
w is —44 for d89l ). well above the height —115 computed above for the leading 
term, see Figure 1. We observe that all of the terms displayed above differ from 
one another either by simple permutations of the wall form components, or by the 
addition of symmetry roots. 

Similar comments apply to the curvature terms with one or more temporal in- 
dices. From d44b . we find, for instance, that terms oc (Co a & c ) 4 or (Coa&c) 2 (Coa<fe) 2 > 
with the indices a, b, c, d, e all different (and no summation on repeated indices!), 
cancel between the two contributions in d44b . Remarkably, these terms are pre- 
cisely of a form which is disallowed by our conjecture above: for instance, if the 
term 

E(R mab f(R mcd ) 2 ~ e~ 2 ° ■ (f+P+r-py (^fi^r-py (99) 

with a, b, c, d all different, had contributed instead of cancelling, it would have 
yielded (choosing indices conveniently) 

w = -(3, 2, 2, 1,1, 1,1, 1,0,0) => C = (0,1,1,2,2,2,2,2,3,3) (100) 

whence £ 2 = +4, in violation of d82l . 

A set of terms which does contribute, by inspection of d44t . is 

&{tt01ab) (-KoiOl) ~ e (101) 

Choosing (a, b) = (9, 10), we get the maximal height contribution 

w = -(1,2,2,2,2,2,2,2,3,3) <=^ w = -[7; 6, 11, 16, 14, 12, 10, 8, 6, 3] 

(102) 

atleveH = —7, of height —93 > —115, and obeying w 2 = —2. The corresponding 
relative vector is 

C = (2,1,1,1,1,1,1,1,0,0) C 2 = 2 (103) 
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and therefore on the boundary of the half-hyperboloid. From the tables of [ 24) we 
deduce that the set of roots corresponding to w belongs to the 5(10 representation 
(010000001) with outer multiplicity /j, = 1. 

There are numerous mixed terms involving the curvature and the 4-form field 
strength, which we again illustrate with some examples (we have checked from 
[43 1 that the structures displayed below actually do appear). For instance, the term 
E(D Fo l23 ) 2 (DoF 0456 ) 2 leads to the £ = 8 root 

w = -(3,3,3,3,2,2,2,2,2,2) w = -[8 ; 5, 10, 15, 12, 10, 8, 6, 4, 2] 

(104) 

with w 2 = —4, and relative vector 

C = (0,0,0,0,1,1,1,1,1,1) ( 2 = 2 (105) 

At level £ = —8, it corresponds to the stio representation (0001000000) with 
[i = 2. The highest root in the multiplet is 

w = -(2,2,2,2,2,2,3,3,3,3) w = -[8 ; 6, 12, 18, 16, 14, 12, 9, 6, 3] 

(106) 

of height —106 > —115, well above the height of the leading R 4 singlet. 

From the C(DF) 3 terms given in 1431 we read off the leading term which is 

ED F 3 4 56 D F 78910 D a F a012 Co bob . The associated root is 

w = -(2, 2, 2, 2, 2, 2, 2, 2, 1, 1) = -[6 ; 4, 8, 12, 10, 8, 6, 4, 2, 1] (107) 

with w 2 = —2, and 

C = (1,1, 1,1, 1,1, 1,1, 2, 2) => C 2 = (108) 

The relevant slio representation at level i = 6 is (000000010) with /x = 1. The 
maximal height in the representation is again above —115. 

Our final example is the combination -E(-Do-Foi23) 2 (Co a oa) 2 - It is associated 
to the £ = 9 root 

w = -(3, 3, 3, 3, 3, 3, 3, 2, 2, 2) = -[9 ; 6, 12, 18, 15, 12, 9, 6, 4, 2] (109) 

with relative vector 

C = (0,0, 0,0, 0,0, 0,1, 1,1) => C 2 = 2 (110) 

Indeed, in agreement with the result d63l above, we see that ( is just an electric 
root. On the other hand, w 2 = —6, and w corresponds to the representation 
(000000100) with fi = 4. The highest root in the multiplet is 

w = -(2,2,2,3,3,3,3,3,3,3) a = -[9 ; 7, 14, 21, 18, 15, 12, 9, 6, 3] 

(111) 

of height —114 > —115, barely above the height of the leading singlet C 4 contri- 
bution. 
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9 Discussion 



Our analysis of higher order corrections to M Theory provides further evidence 
for the validity of the conjecture [ 8 ] that the classical (bosonic) dynamics of M 
theory is 'dual' to a one-dimensional cr-model on the infinite dimensional coset 
space E\q/ K{E\q). In particular, we find it remarkable that the leading wall form 
associated with the known R A corrections, namely the permutation singlet = 
—3(7 does match with a root of E\$ whose corresponding generator is a stio singlet. 
This compatibility between a gravity structure and a Kac-Moody algebra one was 
not a priori guaranteed, and can be viewed as a deep confirmation of the hidden 
role of E\q in M theory. Indeed, as a foil, let us consider the simple case of pure 
gravity in any spacetime dimension D = d + 1. The study of the corresponding 
cosmological billiard has found that one should associate to pure gravity the Kac- 
Moody algebra AEd |4j. Now, for pure (bosonic) gravity, one generally expects 
that the first higher-order corrections will be ~ (R^up^) 2 - However, by using Eq. 
d69b such terms quadratic in curvature correspond to the wall form W2 = — o~ 
whose squared length is (W2) 2 = —d/(d — 1). As the latter squared length is 
never an integer 17 ), we conclude that W2 never corresponds to a root of AEd- 
Let us then consider the problem of determining which values, if any, of the non- 
linearity order N, for corrections ~ R N , might be compatible with the algebraic 
structure of AE&. For instance, let us consider the case d = 3 corresponding to the 
usual (3 + 1) -dimensional Einstein gravity. By using again Eq. d69h one gets the 
wall form wn = — (N — 1) a whose squared length in AEd is — ^zj (N — l) 2 ; 
i.e. — I (N — l) 2 when d = 3. One then concludes that one needs N = 2k + 1 = 
3,5,7,... The lowest candidate for compatibility with AE% is then ~ R 3 . The 
corresponding wall form is W3 = — 2a = — 2 + /3 2 + /? 3 ). Its squared length 
is (W3) 2 = — 6, and W3 is easily seen to correspond to a root a of AE3 of level 
£ = 3 w.r.t. the SL(3) subalgebra [5 ]. However, as in the case discussed above for 
Eiq, the wall form W3 = —2 (/3 1 + (3 2 + /? 3 ) is invariant under permutations of 
the spatial indices. Therefore, for the conjectured correspondence between (3 + 1)- 
dimensional Einstein gravity and AE% to hold, the wall form u>3 should correspond 
to a root of AE3 which parametrizes a singlet of SL(3). However, by using the 
results given in Eq. (8.30) of [5 1 one sees that there is no singlet representation of 
AE3 at level £ = 3. This negative result exemplifies that the compatibility found 
above between ~ R 4 corrections and the presence of 'singlet roots' of Eiq is rather 
non trivial, and could well have failed to hold. 

Our work has also provided evidence for the 'no bounce' behavior of big 

17 For d > 3. Actually, we should restrict ourselves to d > 4 as all R 2 terms are known to be 
on-shell trivial in d = 3 because of the topological nature of the R 2 Euler-Gauss-Bonnet density. 
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crunches in M theory, as naively expected from the dual dynamics, i.e. the global 
structure of null geodesies on E\q/ K(E\q). If we admit the validity of this con- 
jecture, what conclusions can we draw for big crunches in M theory? The 'dual' 
description is an infinite affine length null geodesic going towards f3 — > oo. This 
suggests that the quantization of the E\q/K(E\q) model will exhibit no informa- 
tion loss at the big crunch. On the other hand, the dictionary of [8 9 1 between 
the supergravity description and the coset description is defined only in the quasi- 
classical regime T > Tp. BY contrast, the infinite future of the null geodesic 
motion corresponds to the regime T <^Tp. The fact that the dictionary between 
the two descriptions becomes ill-defined in this limit (and that we find no evidence 
for a bounce) suggests that the infinite 'affine life' near the singularity can only 
be described in the coset variables. This situation is somewhat reminiscent of the 
recent results of [50] based on an AdS/CFT analysis of certain cosmological sin- 
gularities in anti-deSitter solutions of supergravity. 18 [For a contrasting suggestion 
based, similarly to the E\q coset model, on geodesic motion in auxiliary Lorentz 
spaces, see refs. fT6lll7l l. 

However, we feel that at this stage of development of the E\q/ K{E\q) conjec- 
ture, speculating on the ultimate quantum fate of big crunches is premature. Many 
technically challenging tasks remain before one can seriously consider the eventual 
physical consequences of the dual E\q picture. First, one should extend the dic- 
tionary to prove the heretofore unseen roots between height 29 and height 1 15 do 
match in the two descriptions. Second, it would be interesting to explore in more 
detail the validity of the conjecture d82b . which was 'experimentally' observed to 
hold in quite a few specific cases. If our conjecture could be fully verified for the 
known 8-th order corrections, this would open the tantalizing possibility that it still 
holds for higher corrections ~ (R + DF + F 2 ) 7 , {R + DF + F 2 ) 10 , .... It would 
then provide a strong constraint on the algebraic structure of these corrections. As 
these corrections are very difficult to obtain by conventional methods, the hidden 
E\q structure might be of great help in pinning down their structure. 

Finally, among other pressing issues let us also mention: the role of fermions, 
the consequences of compactifying eleven-dimensional spacetime (which is ex- 
pected to reduce the continuous symmetry Eio(R) to the discrete symmetry £io(Z)), 
and the effect of quantizing the coset model (see, in this respect, also the remarks 
in the Appendix below). It would furthermore be interesting to explore the link, if 
any, between the non-compact T-duality symmetries of string theories in presence 
of Killing vectors (which have been shown to survive the addition of higher-order 



The fact that the dual Eio coset picture emerges in full only in the 'strongly coupled limit' 
T <C Tp, while the gravity picture corresponds to the 'weakly coupled limit' T ^> Tp is also 
reminiscent of what happens in the AdS/CFT duality. 
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terms [51 1) and the conjectured Eiq symmetry. 
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Appendix 



Geodesic deviation and sectional curvature on E W /K(E W ) 

As shown in [ 8 , 5 ], geodesic motion on the coset space E\q /K(E\q) is formally in- 
tegrable in the sense that one can exhibit as many constants of motion as degrees of 
freedom (i.e. a doubly infinite tower J) and, moreover, that one can formally write 
a generic solution of the geodesic equation in terms of the constants of motion J. 
We wish, however, to emphasize here that, even if formally integrable, geodesic 
motion on E\q/ K(E\q) is at the same time strongly chaotic in the sense that it 
is extremely sensitive to small deviations in the initial conditions. Let us, indeed, 
consider the 'geodesic deviation equation' which governs the deviation £ between 
two neighbouring geodesies. For simplicity, we focus on the neighbourhood of a 
base geodesic lying entirely in the CSA (corresponding to a simple Kasner solu- 
tion), i.e. where only f3 a and 7r a are excited. This solution is simply described by 
pa _ v aj. _j_ pa w h ere f j s an a ffi n e parameter, (3$ some constants and where the 
'velocity' v a is null: G a f ) v a v b = 0. The other ('off diagonal') variables, u,p, are 
all zero. The geodesic deviation equation governing the affine-parameter evolution 
of the deviation vector £ connecting this geodesic to a neighbouring one can be ob- 
tained in two different ways: either by 'moving' the base geodesic j3 a = v a t + (3$ 
by a generic element in the Lie algebra of Eiq, or by computing the appropriate 
components of the sectional curvature of the symmetric space Eiq /K (Eiq) and 
by writing the general geodesic-deviation (Jacobi) equation. We have checked that 
both methods give the same result. If we decompose the deviation vector £ in its 
various components : 

t = ?ha + 52 tf s) (E^+F^)/V2, 

* — 'a,s V "J 

where h a , a = 1, . . . , 10, is a basis of the CSA, and where £?U denote the compo- 
nents along the non-CSA coset directions. The Jacobi equation yields £ a = 0, 
corresponding to a vanishing sectional curvature in the vv' two-planes (where 
v = v a h a is the vector tangent to the base geodesic and v' = v' a h a a generic 
direction in the CSA), and ^ = + (a(v)) 2 corresponding to a sectional cur- 
vature 

K(v,E+,v,E+) = -(a(v)) 2 

where E+ = (E a + F a )/y/2 is a unit-norm vector in the 'transverse' coset di- 
rection associated to the root a. [For simplicity, we suppressed the degeneracy 
index s on E^ and Fa\] Here a(v) = a a v a denotes the value of the linear 
form a acting on the vector v = v a h a tangent to the geodesic. The point we wish 
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to emphasize here is two-fold: (i) the sectional curvature lZ(v, E+,v, E+) 19 is 
negative (as for geodesies on hyperbolic space), and (ii) this sectional curvature 
decreases without limit towards — oo as the height of the root increases to infinity. 
This result shows that small deviations from the geodesic exponentially increase 
with affine length in a manner which becomes faster and faster as we consider 
components of the deviation in transverse directions corresponding to higher and 
higher roots. This indicates that geodesies on E\q/ K(E\q) are more and more 
exponentially sensitive to initial conditions as one examines them in directions of 
increasing heights. This property might play an important role for the quantization 
of the geodesic motion on Eiq/ 'K(Eiq). Indeed, we recall that quantum motion 
for a scalar particle on curved space involve an ambiguous coupling to the scalar 
curvature, which appears to be given by a formally infinite series for Eiq /K (Eiq). 
This suggests that it is important to include fermions in the Eiq coset model before 
drawing conclusions about the strong sensitivity of classical geodesic motion on 
initial conditions. 



19 As Eio/K(Eio) is a homogeneous space the sectional curvature components computed here 
should be the same in the neighbourhood of an arbitrary geodesic. 
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